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Our universe experienced the accelerated expansion at least twice; an extreme inflationary accel- 
eration in the early universe and the recent mild acceleration. By introducing the Bose-Einstein 
' condensation (BEG) phase of a boson field, we have been developing a unified model of dark energy 

\ (DE) and dark matter (DM) for the later mild acceleration. In this scenario, two phases of BEG 

, (=DE) and normal gas (=DM) transform with each other through BEG phase transition. This 

' unified model has successfully explained the mild acceleration as an attractor. We extend this BEG 

cosmology to the early universe without introducing new ingredients. In this scenario, the infiation 
is naturally initiated by the condensation of the bosons in the huge vacuum energy. This infiation 
and even the cosmic expansion eventually terminates exactly at zero energy density. We call this 
stage as stagfiation. At this stagfiation era, particle production and the decay of BEG take place. 
The former makes the universe turn into the standard hot big bang stage and the latter makes the 
cosmological constant vanishingly small after the infiation. Furthermore, we calculate the density 
fluctuations produced in this model, which turns out to be in the range allowed by the present obser- 
vational data. We also show that the stagflation is quite robust and easily appears when one allows 
negative region of the potential. Further, we comment on the possibility that BEG generation/decay 
series might have continued all the time in the cosmic history from the inflation to present. 

Q_|, PACS numbers: 98.80.-k,98.80.Cq,95.36.+x,95.35.+d 
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I. INTRODUCTION 



Accelerated expansion of the universe seems not to be exceptional. Theoretical and observational studies have 
revealed that the universe experienced the era of accelerated cosmic expansion at least twice; an extreme inflationary 
acceleration in the early universe (EA)[l| and a late mild acceleration (LA)[2] Q. We would like to figure out the 
■ basic physics behind these accelerations and reveal the inevitability, if any, of these accelerations. In the context of 
, the Einstein equation, we need special matter which is endowed with strong negative pressure in order to guarantee 
' the accelerated cosmic expansion. 

\ However, it is often very difficult to find such matter in baryons, and therefore new fields such as the classical scalar 
iy~j Reld has been hypothesized though we have no fundamental understanding of it. Probably reflecting this uncertainty, 
_ this field has been named differently in various contexts in the literature such as Higgs field, inflaton in EA, and 
. K-essence, tachyon, quintessence, Chaplygin gas, phantom field, ghosts, etc. in LA. 

We do not want to be involved in these particular contexts but instead, we would like to reconsider how the degrees 
of freedom of such scalar field can arise within a firm physics even if we still cannot fully identify this field. To this 
aim, we reconfirm the fact that the basic physics which describes the universe from the fundamental level is no doubt 
^ the quantum field theory. Then how the classical dynamical degrees of freedom of scalar field can arise within this 
theory? The most natural mechanism we can consider would be the Bose-Einstein condensation (BEC) of some boson 
field [J] . This condensation process is a phase transition and a typical mechanism that a classical degrees of freedom 
emerges as an order parameter within the quantum field theory . 

Along this line of thought, we have developed a cosmological scenario based on BEC in order to explain LA in the 
contemporary universe^ 0] ■ In this scenario, the dark energy (DE) is identified with the BEC of the bosons 
and the dark matter (DM) with the excited gas of the bosons. This cosmological situation is analogous to the liquid 
'^He system in the laboratory in the sense that the two distinguished components, super and normal, coexist and 
interchange with each other. Our unified model of DE and DM has been successful to describe the contemporary 
universe so far. 

In this paper, we extend this BEC cosmology to the early universe in order to examine how our BEC cosmology 
reveal another acceleration EA. Since we treat LA and EA on an equal footing, we should not introduce any new 
ingredients into the previously developed BEC cosmology except for the cosmological background. Thus the goal 
of this paper is to examine to what extent this BEC cosmology can be verified when applied to the early universe. 
We are especially interested in how the mechanisms of LA and EA are different with each other in BEC cosmology. 
We remember that LA, described in BEC cosmology, has been an attractor Q [1] which guarantees its autonomous 
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realization. Can we also expect EA be an autonomous mechanism in BEC cosmology? 

More precisely, when we analyze EA, we should clarify the autonomous initiation and termination processes of 
the inflationary era, which naturally turns into the reheated radiation dominated era. Actually, the BEC makes the 
order parameter develop from zero to some finite value, and therefore the inflationary dynamics with finite duration 
would initiate autonomously if sufhcient vacuum energy is provided by quantum fluctuations or something. We would 
verify this mechanism first. That is, we derive BEC production rate from the microscopic Lagrangian by using the 
generalized effective action in the in-in formalism. As for the termination of the inflationary era, it may seem to 
be impossible to control the ever increasing order parameter since our BEC model assumes the negative quartic 
interaction (attractive force) without bound. However it is often seen in many physical phenomena that a violently 
unstable mode is taken over by an another stable mode via complex nonlinear processes. We will show that this 
is the case for BEC cosmology. Actually an extreme development of the order parameter completely prevents the 
violent cosmic expansion. We call this process as stagRation, where the order parameter ever inflates while the cosmic 
expansion becomes stagnant. Thus the inflationary era terminates autonomously within a finite time scale. Moreover, 
this extreme development of the ever accelerating order parameter easily gives rise to the particle production, which 
provides a natural process of the reheating. Furthermore, the instability analysis reveals that the homogeneous mode 
immediately decays into the inhomogeneous localized modes at this stagflation era, when the cosmic expansion ceases. 

This stagflation era also provides a key mechanism to solve the cosmological constant problem. The issue is to 
explain the dynamics of the cosmological 'constant' A, which had once been a huge value (10^^~^^ eV)^ that makes 
EA, turns into a tiny but non-zero value about (10~^^ eV)^ that makes LA. Though this initial huge value for A 
would be naturally provided from ubiquitous vacuum energy, there should be a certain mechanism for the very fine 
tuning of A to the vanishingly small but finite value at present. In the previous analysis of BEC cosmology applied 
to LA, we assumed that the zero-point of the potential energy is exactly zero, and shown that the vanishingly small 
but finite value of the effective A term appears as a fixed point, in which the condensation speed and the potential 
force balance with each other. In this paper, we propose a new mechanism which autonomously adjust the effective 
A term exactly zero. This would justify the above assumption for LA. The essence of this A-extinction mechanism is 
the fact that the stagflation takes place exactly at the zero total energy. As explained above, the uniform mode of 
BEC becomes unstable at this stagflation point and suddenly decays into localized objects. Thus after the inflation, 
the universe evolves into the stagflation point, where A becomes zero and then the universe becomes dominated by 
the reheated radiation. We would like to quantitatively verify this mechanism in this paper. It will turn out that the 
A-extinction is asymptotically realized after many BEC decays. 

The above instability of the uniform mode is the same as that appears in BEC cosmology application in LA. 
Moreover, the condensation process of BEC commonly triggers all the essential dynamics of the cosmic history in 
both the cases of EA and LA. It is important to notice that these decay due to instability and the condensation of 
BEC are not a coherent time change described by any Lagrangian but are rather phase transitions and the incoherent 
evolution with dissipation and fluctuations. This point is further discussed in the following sections. 

This paper is organized as follows. After applying the in-in formalism into BEC formation, we review the essence 
of BEC cosmology applied in LA in section 2. In section 3, we apply the BEC cosmology to the early universe in 
the order of initiation, termination of inflation, reheating, cosmological constant problem, and generation of density 
fluctuations. In the last section 4, we summarize our study emphasizing the robust stagflation and clarifying unsolved 
problems at present. 

II. BASICS OF THE BEC COSMOLOGY 
A. BEC in the Universe 

In this section, we first review our model of BEC cosmology from the microscopic physics. Various observations and 
theoretical analysis have revealed that the DE and DM dominate the matter contents of the universe and DE/DM 
ratio is of order unity at present @. These facts may indicate some deep relation between DE and DM. Therefore it 
seems natural to construct a unified model of DE and DM based on a single physical principle. The simplest thought 
will be to consider that both DE and DM are actually the same dynamical matter but in different phases of existence. 
DE phase should be almost uniformly distributed and is supposed to provide accelerated cosmic expansion. On the 
other hand, DM phase should be inhomogeneously distributed and is supposed to trace the structure of galaxies and 
clusters. 

According to the Einstein equation, the acceleration of the cosmic expansion requires isotropic negative pressure 
comparable to the energy density p < -~p/3 for DE. This is because, in general relativity, pressure as well as energy 
density works as a source of gravity. Such strong negative pressure sounds strange if we consider fermionic fields which 

constitute ordinary baryonic matter. On the other hand for bosons, familiar static magnetic field H , for example. 
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shows such behavior p = —p = —H'^/ (Stt) in one spatial direction, reflecting that the magnetic force hne has a 
positive tension. This negative pressure would be the origin of the variety of bursts and accelerations in solar physics 
such as flares and coronas. Unfortunately this is a vector fleld and shows positive pressure in the other two spatial 
directions. A natural possibility to obtain the isotropic negative pressure is to consider the classical scalar field. This 
would be the standard line of thought in the modern cosmology as explained in the introduction. 

However, when we construct a unified model of DE and DM, we further have to specify the origin of such scalar 
field in order to clarify the connection between DE and DM. The problem is to find a phase which behaves as classical 
scalar field in the matter which is originally described by quantum physics. The most natural case would be the 
quantum condensation of bosons in low temperature, i.e. the Bose- Einstein condensation. Macroscopic number of 
boson particles share a single ground state in common and behaves as a coherent wave. This condensed phase is well 
described by a classical mean field. On the other hand the rest of the boson gas in excited states should macroscopically 
behaves as classical gas. Therefore we naturally arrive at the thought that the boson gas is composed of two phases, 
condensed phase and gas phase, each of which can be regarded as DE and DM, respectively. Thus the BEG cosmology 
arises. Since the basic structure of the universe in this model is quantum physical and cosmological, we have a unique 

mass scale m — {^H^h^/ (^Gc^))^^^ = 0.0092eV obtained by equating the present energy density p = H^/G to the 
energy density associated with mass m and the Compton wave length I — h/ (mc). This scale characterizes the cosmic 
BEG. 

The above two-phase structure is quite common in physics in the sense that any special interactions among boson 
particles are not required for BEG and it inevitably takes place provided the temperature is lower than the critical 
one. A well known such structure would be the liquid ^iJe system in the laboratory. The super phase and normal 
phase correspond, respectively, to DE and DM; they coexist and interchange with each other. Although the origin of 
the two-phase is quantum phase transition, the condensed phase is a classical mean field and should not be confused 
with the macroscopic quantum state, which is described by a wave function. 

The BEG of boson mass m takes place provided the de Broglie wave length A^s = 2TTh? / [mkT) exceeds the mean 
separation n^/^ 



kT < kTc; 



(1) 



This condition requires the boson gas to be 'cold', if we were to realize BEG for reasonable boson mass appropriate 
for DM. Actually in the universe, if the boson is the cold dark matter, irrespective of its interaction, its temperature 
(or quasi-temperature) behaves as the same as above T cx n^^'^ in the cosmic expansion. Therefore BEG can always 
take place and continue provided the boson temperature was less than the critical temperature at some moment of 
cosmic evolution. The corresponding critical boson mass for BEG can be arbitrary in general, but if we temporary 
assume that the boson was in equilibrium with the rest of the matter in the past, then the upper limit of the mass 
turns out to be about 19eV [8]. 

In the universe the ratio of DE and DM cannot simply be given by the standard expression for the ratio of 
condensation 1 — {T /Tcr)^^"^ in equilibrium with the temperature T. The universe is actually almost adiabatic and 
therefore the transition toward the ground state is forbidden in the first approximation. However a small non- 
adiabaticity allows the very slow condensation. We need a firm formulation which describes this condensation process 
based on microscopic physics. The most appropriate, but not necessarily complete, method for this purpose would 
be the generalized effective action formalism Q, [lo| . [Tlj . on which we now briefly explain. 

The partition function, which describes all the fluctuations, is given by 



Z[J] = tr 



T 



exp 



J(t> 



exp iW[J] 



(2) 



where the integration is taken over the entire space and the doubled time contour which extends from t = —oo toward 
t — +00 and then back to t — — oo again. The symbol T rearranges the operators on its right hand side in the order 
of this doubled time contour. All the quantities with overtilde represents that they are defined on this contour. J is 
the external source field, p is the initial density matrix, and the quantum field cf) is in the Hcisenberg representation. 
This partition function can also be represented in the interaction picture as 
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Z[J] = tr 



T exp 



i Jc^-i V{cj>) 



exp 



tr 







" 6 

_i6J 




exp 


HI 


exp 




^ J0 







J (x) Go {x,y) J (y) 



(3) 
(4) 



where the Lagrangian is decomposed as the free part Lo[(j)] and the interaction part V[(j)], L[(j)] = LQ[(j)] — [</>]. 
The Wick theorem is used for the last equaUty. It is useful to introduce a vector and matrix representation for the 
functions on the generalized time contour; a field X{x) can be written as X^{x) + where the quantity with 

suffix +(— ) has its support only on the forward (respectively, backward) time contour. Then X{x) can be represented 
as a vector For example J (x) = {J+{x), J-{x)) and (j){x) = 0_ (a;)) . We also use the 

notation Xa{x) = X-^-{x) — X^{x) and Xc{x) = (X_|_(.t) + X^{x)) /2 in the following. Then the notation / J(p means 
J J(j) = J dxJ-^-{x)4i^{x) — J dxJ^{x)(l)_{x), where the subtraction on the right hand side represents the fact that 
the direction is opposite in the backward time contour. This can be conveniently expressed in the two-by-two matrix 
notation, 



J(j)= dx(J+{x),J^(x))ai{(l)^{x),(j)_{x)f 



(5) 



1 



where "'3 = ^ q ^ j transpose. 

The generalized two-point function GQ{x,y) = ~i (T(l){x)(j){y)) is also expressed in this matrix notation as 



Gf ix,y) G+ {x,y) 
G- {x,y) Gp {x,y) 



-I {mx)cj){y)) 
-i {<j){x)(f>{y)) 



-i {(j){y)(j){x)) 
-z(r^(x)0(y)) 



(6) 



The fact that the only three of them are independent is most clearly observed in their momentum representations. 



GF{k) 



-D{k)+iB(k) 



Gpik) 



-D{k) -iB{k) 



,G±{k) 



A{k)TB{k) 



(7) 



D(fc)2 -f A(fc)2 ' ' D{kf + A{k)^' ' D{k)^ + A{kf' 

where A{k) ,B (fc) , and D (fc) are mutually independent functions. 

If we were to define the classical field as ifix) = SW/SJ and define the Legendre transform T[ip] = W[J] — J J if 
as the effective action as usual, then by definition, the equality holds 5T /5ip{x) — —J{x), which may be regarded as 
the equation of motion for ip{x) including all the quantum fluctuations. However, reflecting the fact that the effective 
action is not in general real T = ReT + ilmV, the above method makes the field 1^(2;) unphysically complex. This 
ordinary procedure cannot describe the phase transition dynamics. Therefore we go beyond the standard method and 
decompose the effective action. We first notice that ImT is even in the variable (p^{x) = ip^{x) — ip_{x), 



ImT 



ip^{x)B{x - y)(p^{y)dxdy 



(8) 



and the kernel B{x — y) is positive definite[9']. Therefore, introducing an auxiliary real field f (x), we have the following 
path-integral expression. 



= / [dC]P[£,]eTip[iReT + (x) (p ^ (x) dx 



where the weight function is defined as 



i{x)B-\x-y)i{y)dxdy], 



(9) 



(10) 



where N is the normalization factor. Thus the field <^ [x) can be regarded as Gaussian random field. Inclusion of 
higher order terms in the above makes us possible to treat the non-Gaussian properties as well. On the other hand, 
the real part of the action 



ReT 



(11) 
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should be regarded as the standard action which describes dynamics. Applying the variational principle for the 
variable ipi^{x) 



= -Jc{x) (12) 



we obtain the evolution equation for ^Pci^) = {f+i^) + '^^-i^^) /2, which is simply denoted as ^p{x) hereafter, 



{dd + m^)ip+ dt' dx'A{x~x')ip{x')=i{x), (13) 



where the source term Jc{x) is dropped. This is the Langevin equation for the classical order parameter ip{x) with 
classical random field ^ {x). Actually, if we define the statistical average by 

(...), ^ j m...p[s], (14) 

then the correlation function for the auxiliary fields becomes 

{ax)ay))i^B{x-y). (15) 

Moreover, Eq. (fT3|) is necessarily causal 9] since the kernel A{x — x') has the property A{x — x') — for < 
Thus the full effective action is now rewritten as the statistical superposition of dynamical effective action. Note that 
the above dissipative nature arises as a result of an implicit and systematic coarse graining of the original full theory 
through the process of renormalization in the in-in formalism. This renormalization process is generally inevitable in 
quantum field theory to obtain finite physical quantities. The dissipative nature however cannot arise in the ordinary 
vacuum as in the standard quantum field theory. For dissipation to arise, the system must be unstable or permit 
particle production or be immersed in the dissipative environment. 

The above Langevin equation can, in principle, be transformed into the equivalent Fokker-Plank equation for the 
distribution function. This process is actually tractable if we can use the Markovian property. 

The most of the application of the above formalism to the BEC is given in [ll|. The Langevin equation is transformed 
into the Fokker-Plank equation 

for the probability distribution function P[(j);t\, where e' and fi are, respectively, renormalized energy and chemical 
potential which are related to the function D, and R, are the transport coefficient functions which are related 
to the function A and B. This Fokker-Plank equation is further transformed into the rate equation for the energy 
density of the condensation p (fc; t) . In the limit of small condensation, i.e. the back reaction of the condensation to 
the gas is negligible, the rate equation reduces to 

lpik;t)^r-(i-^^)^r-, (17) 



where p^'^ {k;t) is the equilibrium value of the condensation and F™ = —R + iY,^/2. The order estimate of the results 
in [Tli becomes, 



if the dominant interaction is A {(fxj)*)'^ , and 



\\fm'^pl,,, (18) 



r»«HVp^,„ (19) 



if the dominant interaction is of type a(j)Ad4>* , where A is the radiation field. In either case, the induced process 
is neglected in the RHS of Eg. p?)) . It is most interesting that the detail of the form of F*" is not relevant for the 
BEC cosmology (fs^). Any choice of F™ (or F in the later argument) yields qualitatively the same cosmology. In this 
fiexiblc possibilities, we adopt Eq. fTS)) in the following although Eg. p^ is equally applicable. 
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The above results determine the BEC condensation rate in the uniform boson gas. Note that the above condensation 
rates are transport coefficients which cannot be found directly in any Lagrangian. This fact yields essential difference 
between our approach and the presently popular interacting DE-DM scalar field models. For example [l^ 



study DE-DM coupling as an interaction term in Lagrangian. The recent reference [IJ] study DE-DM coupling 
phenomenologically by setting w = p/p as a constant parameter. The quantity w actually evolves in time in our 
model according to the equation of motion. In these aspects, our BEC model should be distinguished from general 
DE-DM coupling models. 

So far we have considered the generation of BEC. While this condensation proceeds, the uniform BEC eventually 
becomes unstable due to the Jeans instability and collapses into localized objects (inhomogeneous DM) leaving no 
uniform component. If the boson gas has attractive interaction, then this collapse would be further enhanced. 

Thus in the universe, BEC is generated and collapses. These dynamics can be regarded as generalized phase 
transitions which cannot be described only by the coherent time change governed by a Lagrangian. We now examine 
the detail of the BEC phase transition in the following subsections. 



B. Homogeneous mode: steady BEC 

The gradual generation and subsequent decay of the spatially uniform BEC can be summarized by the following 
set of equations in a compact form 

^ = [a) -^iP.+P^ + P^)^ 
■p^ - -QH{p^-V)+Tpg-T'p^, 

■pi - -iHpi+V'p^. (20) 

where we assume a simple form of potential V = Vo + m'^(f>*(p + ■|((^*(^)^ with attractive interaction A < 0. This 
potential form will be fully fixed if we can identify the boson field in the future. The Einstein equation in the first line 
of Eg. ((20)) gives the cosmic expansion rate H in terms of the energy densities of gas Pg, BEC p^, and the localized 
objects Pi- The second line represents the evolution of the gas energy density Pg which reduces due to the cosmic 
expansion —SHpg and the consumption for BEC —^Pg- The third line is the relativistic form of Gross-Pitaevski 
equation describing BEC dynamics with the steady condensation rate Lp^, where F = const, (for Eq. (jl9p ) or F cx 
(for Eg. lfTS)) ). This form of source and loss terms ±Fpg is the most characteristic to our BEC model although they 
are small. In the same line of Eq. pp)) . the term — F'p^ symbolically represents a sudden collapse of BEC. The detail 
of this collapse will be explained in the next subsection. The last line represents the evolution of the localized energy 
density pi which reduces due to the cosmic expansion —2>Hpi and increases due to the BEC collapse +T' p^. We now 
explain the BEC instability and the collapse. 



C. Inhomogeneous modes: BEC collapse 

Uniform BEC is not an absolutely stable phase if the interaction is attractive or gravity is concerned. It collapses 
into localized object after the appearance of instability. Let us estimate the scale of BEC instability, which determines 
the time scale of the BEC collapse and the resultant mass scale of the localized objects. We start the space-time 
metric of the form 

ds^ = (1 + 2$)dt2 - 0^(1 - 2$)dx^ (21) 

where $ = $ (i, x) represents the gravitational potential and a = a{t) the cosmic scale factor. BEC is described by 
the equation of motion for the order parameter (p ~ (j){t,x), i.e. the third line of Eq. ipDl) neglecting F and F' since the 
instability time scale is very short compared with F^^ and F'^^. Up to the first order in $, this equation becomes 

ij) + 3-0 - 4tV2(/) + m^(l + 2$)0 + A(l + 2$)|(/)p(/) = . (22) 
a 

Gravity is described by the Poisson equation, in the Einstein equation, 

V^* = 4nGa' + • -0^ + m'\cP\' + ^101' - Pol , (23) 

I a a 2 1 
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where Pq is the uniform background energy density. We consider the linear perturbation on the uniform background 

(l>it,x) = ^Q{t)+^^{t,x) ,^{t,x)^0 + ^i{t,x). (24) 
where small perturbations (f)i and $i are decomposed into Fourier modes, 

(pi — (l)iexp{nt + i k^), 

$1 = $iexp(m + «'fc^). (25) 
The background variables satisfy, from Eqs. (1221123]) . 

00 + 3^00+ K + Al^oH-^o = 0, (26) 
■ ^rn^ + Tl'l^oAl^ol'-Po = 0. (27) 



4' 



In these Eqs. (|26ll27|) . the variable (pQ can be assumed to be real without generality since all the coefficients are real. 
Then the existence of a non-trivial solution requires that the above set of linear equations are dependent with each 
other. Thus we have the condition 

a / 47rGa^ 



2(m^ 



X4>l)cPo (20of^ + 2(m2+A<^g)0o) 



X [fl'' + 3-n + %+m' + \<j)f,] =0, (28) 
\ a a-^ J 

which determines the instability strength 17 as a function of the wave number k. The wave number k which corresponds 
to the maximum value of (i.e. most unstable mode) is not necessarily the relevant scale since we also have to consider 
the completion time of the collapse a/ (ak), where a (< 1) is the decay speed of BEC. Thus the relevant scale will 
be the wave number fc, which should satisfy the condition that the collapse time is just the instability time i.e. 
afc,/a = Q. Choosing the value associated with the most unstable mode among four solutions of Q, we have 

1/2 

(29) 



f.^ / -miff + ^mjff + 647rG(l + a^) (m^^ _ 2Km^4 + K^cj>^„) ' 



a 



2(1 + a2) 



where miff = m^ — 3k0q, k = — A(> 0) and adiabatic approximation _ff = 0, 0o = is utilized since the collapse time 
scale is much smaller than cosmic and condensation time scales in the late universe. 

The actual scale critically depends on the boson massQ. If it is about leV, then the typical mass of the structure will 
be M, « 1.6 X IO^^Mq which is of order of a galaxy. If it is about 10~'^eV, then the typical scale in tiny M* « 10~^g 
and therefore BEC collapse cannot form macroscopic cosmological structures. In any case, these collapsed clumps 
would work as the ordinary DM and the scenario of large scale structure formation exactly reduces to the standard 
model. 

A typical numerical demonstration of uniform BEC mode is given in Fig.([T]). Most of the characteristic features of 
BEC cosmology are shown in this figure. The gas density Pg (grey line) reduces monotonically while the BEC density 
P0 (solid line) gradually increases toward some instability point where it suddenly decays into localized objects pi 
(dashed line). Then the BEC density disappears and a new BEC begins. This cycle repeats several times and finally 
the phase of p^ = const, appears. Thus the universe inevitably enters into the final accelerated expansion phase. 
This phase is guaranteed by the stable balance of the condensation flow and the potential force, and is an attractor 
of Eg. ((20)) . Several other features and predictions of this model are found in [8]. 



III. BEC IN THE EARLY UNIVERSE 



We now extend the BEC cosmology to the early universe. As was studied previously, for the relativistic case, the 
critical temperature T^r oc p^^"^ increases much faster than the cosmic temperature T cx p^/* in the course of increasing 
energy density p. Therefore we safely expect the universe is always under critical temperature provided so at present. 
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FIG. 1: A typical numerical demonstration of the BEG cosmology. This graph represents the evolution of energy densities of 
the boson gas (grey line), BEG (solid line), and the localized objects p; (dashed line) as a function of redshift z. We solved 
Ea. (l20|l with the parameters A = -0.51, V = 0.2601, m = 0.002861^. The BEG GoUpase time sequence is z ^ 69, 54, 39, 26, 15. 
After four collapses of BEG, the energy density of BEG p^ becomes constant, which induces the accelerated expansion of the 
present universe. 



The acceleration in the early universe is quite different from that in late universe in the sense that a huge vacuum 
energy Vb is present from the beginning in the former. However we do not want to introduce any new ingredients into 
the model used in the original analysis in LA. This is because we want to understand the cosmic accelerations from 
a unified point of view. Actually, the mass m of the boson is about eV and therefore is completely negligible in the 
early universe. The condensation rate T has been quite robust in the previous study in the contemporary universe 
and therefore we take the same form of F for BEC condensation as used previously. 

The goal of this section is to examine the following individual steps to complete the whole BEC cosmological scenario 
and verify this model: 1. Natural initiation of the inflation from the fireball stage is expected from the steady BEC 
processes. 2. This inflation stage terminates autonomously due to the stagflation. 3. Reheating process inevitably 
takes place due to the non-linear coupling in BEC. 4. The cosmological constant is autonomously adjusted to be zero 
due to the BEC instability at the stagflation point. 5. We calculate the generated density fluctuations and compare 
them with observations. 



A. Natural initiation of the inflation 



The original model of the BEC cosmology is described by Eq. (|20|) . Now in the application to the early universe 
whose energy scale is extremely higher than eV, a typical mass scale of the boson, relativistic damping terms should 
be used: 

Pg = ~4:Hpg-Tpg, 

P<p = - V') + Vpg - F'p^, 

■pi = -AHpi+T'p^. (30) 
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where = 87rG'/(3c^), the potential is V = Vq + ni?(j)*(j) + with A < 0, and the energy density for BEC is 

= 0*0 + V. 

In the early universe, tiny mass term is neglected. On the other hand, the condensation rate F, whose energy 
scale dependence is not clear, may still be relevant in the early universe. We should consider an initial vacuum energy 
Vq ~ mpj, which is huge up to 10^^° times the present value. 

Initial boson gas Pg (t) = p^j {ti/tf' and the temperature will be reduced by the cosmic expansion a{t) = oq {t/ti)^^'^ . 
On this background, BEC gradually proceeds since Tcr < T is guaranteed from the energy density dependence of 
them as argued above. Since the initial evolution velocity of the BEC (f) is small and the potential energy is negligible 
compared with the boson gas density, the BEC approximately obeys 

3i70'=rpg, (31) 
whose solution, under the choice F (i) = ^oPg {t), is given by the fire ball solution. 

The situation of each stage of BEC will be given in Fig. ([4]). Then, if we identify the initiation time of inflation tunf 
by the condition 

Pg{tnn{)^Vo, (33) 

then the initial value of the BEC when the inflation begins will be 



3/i 

or 



H.(W) = 2t.(£^#!V'7l-f^V'1 (34) 

9« 



^-0..e(iiin.)=f^^V'V^V''- (35) 



3/i 



Although the boson gas will be diluted, this initial BEC triggers the further evolution of BEC, and BEC proceeds 
consuming the vacuum energy. During the inflation, BEC evolves as 



which solves, under the initial condition Eq. (j35p . as 

1/2 / 



-2A 



iiinf — t + 



1 fm. 



2 \-XS' 



-1/2 



(36) 



(37) 



where we approximate = ^inf = M^^o- Since the potential is flat around the origin, this solution actually satisfles 
the slow-role condition 











«1, 


^inf 


2V 


3i?inf0inf 



< 1 



or equivalently under the condition Eq. p6p . 



|£| < 1, |?7| < 1, 



where 



ml, (V^ 2 



167r V V 



77 ' ^ 



(38) 



(39) 



(40) 
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These conditions hold until just before 0inf (0 tii^s the (virtual) singularity a,t t = toinf (where (j)ini diverges). Thus 
the duration of the inflation is 

1/3 

Ainf ^ teinf " ^iinf = ^ + O (ffr^) (41) 

and the e-folding number N becomes N — i/inf Ainf. 

In Fig. ([2]), we have depicted a typical time evolution of the scale factor. We do not seek for the most favorable 
values of parameters in this paper since our aim is to clarify the basic mechanism and to show the robustness of the 
BEC cosmological model in the early universe. 




5 10 15 20 25 30 35 

t 

FIG. 2: The time evolution of the scale factor a (t) . The scale factor is normalized by itself at the begining of the inflation. 
This is a solution of Eq.(l30} with parameters: A = -0.5, Vb = 1, 4>{0) = 0, </> (0) = 0.001, .pg (0) = 5, To = 1, in the units 
of Hini for H, H[^f for t, Hi^-^i/fi'^ for p, where = p^Vb. The inflationary phase, i.e. a oc e^'"'*, naturally appears with 
sufficient e-folding number. The total energy density vanishes and the expansion ceases at the stagflation point t — ts ^ 34.8, 
where the strong instability of the uniform BEC mode actually invalidates the evolution after ta . The inset magnifies the local 
region around ts, which is marked by an arrow. 

This mechanism of inflation is the same as the standard models though the initiation of it is quite different. 
Moreover, the termination of it is also different as we will see in the next subsection. In the inflationary phase, BEC 
develops by consuming the vacuum energy contrary to the case of BEC evolution in LA in which BEC develops by 
consuming the boson gas energy. 

B. Termination of inflation 

As is shown in Fig. ([2]), this inflation sharply terminates within a finite time. This can be easily understood as a 
general feature of our BEC cosmology. In the late stage of inflation F, F', are neglected, and Eg. ([30]) reduces to 

H' = ^P,,P, = -3H^\ (42) 

and, provided H 0, further to 



(43) 
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Since the right hand side of Eq. ([43|) is negative and (p increases in time, the expansion rate H eventually crosses 0, at 
time t = ts, after that the universe terns into the contraction phase. The dynamics of BEC in this process corresponds 
to the over-hill regime in the original BEC cosmology [7]. 

When the cosmic expansion ceases at t = tg, Eq. pO]) yields, neglecting Pg, 

H = Q,p^ = Q,p^^Q, (44) 

which guarantee the energy density smoothly vanishes at t = t^. However, we should notice that H from Eg. ([^5]) 
and the continuity. We call this phase as stagflation regime, i.e. the cosmic expansion is 'stagnant' while the BEC 
is steadily accelerating and 'inflating'. In the previous approximate solution Eq. (|37p . the termination time tdnf is 

almost tg-. ts toinf- 

Figure ([3]) shows the energy densities of the boson gas Pg (grey line), and of BEC p^ (solid line) in the same time 
evolution as in Fig. ([2]). The stagflation regime is magnified in the inset of this figure. The dashed line shows the 
Hubble parameter which smoothly changes its signature, expanding to contracting at tg. 




FIG. 3: This graph represents the evolution of energy densities of the boson gas (grey line), and of BEC p^ (solid line) as 
a function of time t. The enrgy density is normalized by {Hint/p)^ and the time by Hr^j. The inflationary phase eventually 
terminates and the universe turns into the contraction phase. It finally approaches the virtual singularity, if we do not consider 
the instability and the decay of BEC. Inset is the magnified graph at the last stage of the infiation. Dashed line is the Hubble 
parameter, which crosses zero and the total energy density vanishes at the stagflation point t — ts ~ 34.74 (marked by an 
arrow). The strong instability of the uniform BEC mode actually invalidates the evolution after ts- 



We would like to comment briefly on the virtual history of the universe after the stagflation regime though it 
actually does not exist, as will be explained in the following subsections. The stagflation regime Eq. ((44)l docs not 

• 2 4 

last forever; the kinetic term (f> eventually develops and surpasses Xif) term. In this kinetic term-dominant regime, 
Eq.® yields d (^0^2) /dt = -?>Hip^ and iJ^ = p?ip^ 12, and therefore 

3-2 

(/) = -j^pc^ , (45) 
where we must choose H < and (f> > for consistent solution. These lead to the implosive solution 

3p 

where the universe hits the singularity at time t^. The cosmic contraction becomes 

a = const, (i, - t)^/^ , (47) 



imp = const. - — log (i, - t) , (46) 
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which explains the virtual behavior in Fig. ([2]) for t t^,.{See also Fig. ([4]).) 

As a whole, the early universe suffers the BEC phase transition as depicted in Fig. ([3]). This figure represents the 
whole evolution of BEC field numerically obtained, and several analytical approximations whose validity is 

limited within each local region. Analytic solutions are 0^^.^, in Eg. ([5^ . cji^^f in Eq. ([57|) . ^gt^g Eq. ([5T|l (which will 
be described shortly), and 0i,„p in Eq. (j46[) . The solutions ^gj.^ and (p^^f are matched at tiinf, defined by Eq. (j33p . The 
stagflation time ts, marked by an arrow, is defined by Eq. (l44|) . The implosion time is the time when BEC field 
diverges as Eq. (|i5|) . 




FIG. 4: This figure represents the evolution of BEC field 4>{t), numerically obtained, and several analytical approximations, 
'/'firo) 4'infy '/^stagi ^-^d (jtimp- The applicability range for each solution is shown by shaded region. The solutions <j>firay 4'inf ^-^e 
matched together at tuni ~ 0.16. The stagflation time ts ~ 34.7 is marked by an arrow. The BEC field (j) diverges at the 
implosion time t* « 35.1. 

In the real world however, the infinite amplitude (f) oo and velocity (p —^ oo for BEC are not allowed. Actually, 
before the universe enters into this phase, there appear two new features in this scenario and the singularity, which we 
call implosion, is avoided. One new feature is the particle production due to the ever accelerating BEC (p oo. This 
is the origin of the reheating that some portion of the original vacuum energy Vq is converted into thermalized gas 
of bosons. Another new feature is the instability of the uniform mode of BEC at t « and ~ 0, and subsequent 
decay of the uniform BEC into the localized objects, without leaving any vacuum energy. This is the mechanism 
which did adjust the final vanishing vacuum energy, or vanishing cosmological 'constant', despite the huge vacuum 
energy set before inflation. Let us examine each of these features separately in the following. 



C. Initiation of reheating 

The rapid increase of BEC ^ oo at the late stage of the inflation yields plenty of particle production, which will 
reheat the universe, through the self-interaction. If this particle production is effective to turn sufhcient amount of 
kinetic energy of BEC into thermal gas, the universe can go into the big bang phase. In order to demonstrate this 

• 2 -2 

scenario, we consider a simple model of particle production; an addition of the term —"ftp (+7^ ) to the right hand 
side of p0 {pg, respectively) in Eq. pO]) . We can approximately evaluate the particle production during the inflationary 
era. By using Eg. ([57)1 . we obtain 
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If we set 7 « w leV, then the reheating temperature, using the expression ~ ^ ) ^o^l'i be 

/lO \^^^ 

ksTrh < —imli (Hinf/m^)'/^ « 1.17 x lO^^ (Fi„f /mp/)'/^ GeV. (49) 



Note that the particle production in the BEC model is caused simply by the accelerated evolution of the background 
field <l). 

D. Autonomous adjustment toward the vanishing cosmological constant 

Let us discuss on another new feature which does not allow the implosive singularity. As in the previous arguments 
on the BEC cosmology in LA, the instability of BEC actually develops and the homogeneous BEC phase collapses 
into many localized objects. Thus dynamics of BEC is not closed within the homogeneous mode contrary to the 
recent arguments on the DE-DM coupling ([l3|) for example. 

In the middle of the inflationary phase, the BEC fluctuations cannot develop. However, they have a chance to grow 
after the inflationary phase. We now examine this around the stagflation time when the energy density vanishes at 

t = ts. 

The middle of Eg. ipi)) allows us to approximate the evolution equation, neglecting V07 as 

<?^'=^/^0^ (50) 

which immediately solves 

-^stag = {ts - , (51) 

where k = — A(> 0). In this regime t ^ tg, the instability scale Eg. ipS)) becomes 

^* fa ^l/4/-<l/4 1/2,3/2 /„ -1 (. +\-3/2 /rrjN 

— = (Stt) ' G ' k ' 03t'ag = (Stt) ' G ' k (ts-t) ' , (52) 

which increases without bound for t tg. On the other hand, the cosmic expansion is almost ceased in this regime 
t Ki tg. Therefore, we expect a strong instability takes place and BEC rapidly collapses into localized objects, whose 
density is denoted as pi- Moreover, there is sufficient time {tg — t) for the completion of the BEC collapse before the 
universe escapes from this regime t Ki tg-. 

{tg ~t)> ajK oc {tg - tfl"^ . (53) 

This is similar to what happens in the over-hill regime in the late time BEC cosmology. Thus all the BEC component 
turns into localized objects, which behave as dust fluid, leaving vanishing vacuum energy or cosmological 'constant', 
despite the initial value of the vacuum energy Vq. The fate of these localized objects is not yet clear. If they form 
light blackholes, they may eventually evaporate into radiation including the boson gas. If they form heavy blackholes, 
they may remain until now and yield significant effects on the cosmic structure. We simply choose the former case in 
our numerical demonstration in this paper. 

More precisely, we have numerically checked the above instability directly based on Eq. (f28|) without the approx- 
imation Eq. (|^5)) . It has turned out that the collapse of BEC, estimated from the linear instability analysis, takes 
place just before the stagflation point tg. This means that a small amount of vacuum energy remains at this point 
since p = realizes only exactly at i = is. However further evolution of the uniform BEC from that point leads 
to the next instability of this BEC, which yields much smaller remnant vacuum energy. These sequence of BEC 
collapse continues many times and the vacuum energy reduces at each time. Figs.([5][Zl) show a typical example of 
the numerical demonstration of this sequence. Fig. ([5]) shows the time evolution of the BEC order parameter (/>. Each 
segment represents a history from the generation toward the decay of BEC. The average value of reduces in time 
and the duration of each history becomes longer. Fig.® shows the time evolution of the scale factor a(t). Initial 
inflationary era, which is a straight line on this graph, is followed by much milder expansion era. The latter is shown 
in the inset of this figure and is almost the same as radiation dominated cosmic expansion. Fig. ([7]) shows the time 
evolution of BEC energy density (solid line), boson gas energy density p^ (gray line), and the Hubble parameter H 
(dashed line). The reduction sequence is apparent. In our demonstration Fig.©, the vacuum energy at each segment 
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era reduces as 1.0, 0.0063, 0.00072, 0.00016, 0.000046, 0.000014, 4.6 x 10"^ 1.5 x 10"^... in the unit of the original 
vacuum energy Vq. 

The above special feature mainly comes from the stagflation, which induce the decay of the uniform mode of BEC 
or classical scalar field. The stagflation itself can be realized when the potential allows negative region, which will 
be further discussed in the summary section. Thus the stagflation and the associated reduction of the cosmological 
constant are not special features in the universe. 



1.50 
1.00 
0.70 
0.50 




FIG. 5: Time evolution of the BEC order parameter (j> is shown. It shows a sequence of BEC collapses. It may happen that 
this BEC reduction sequence leads to the era of LA in Ref.[T]. However, the numerical calculation extending over 100 digits is 
almost impossible and we could not reveal the full extention of the reduction sequence. 



E. Generation of density fluctuations - comparison with observations 

As we have examined in the previous two subsections, both the reheating due to the particle production and the 
decay of the vacuum energy due to the BEC instability are most effective at around the stagflation time t = ts- 
Furthermore, we observed that this set of BEC generation and decay repeats many times and this sequence makes 
the effective cosmological constant reduce toward zero. Although the detail depends on the BEC decay dynamics 
which is quite nonlinear and goes beyond the present paper, we may naturally assume that the universe approaches 
toward the hot FRW universe. Then its fate will be well described by the standard cosmology until it enters into LA 
regime in most recent era. Therefore the most characteristic feature of our model will be related to the phenomena 
around the stagflation point tg . Then one of the remnant characteristic features to our BEC model will be the power 
spectrum of the produced density fluctuations reflecting ever increasing development of the BEC component toward 
the end of the inflationary era. 

We now calculate the produced density fluctuations during the inflation. Since our BEC behaves almost like scalar 
field during the inflation and therefore we can safely apply the standard analysis on the density fluctuations. Based 
on the squeezed quantum state in de Sitter space with slow-rolling approximation, we can estimate the curvature 
fluctuations and the tensor fluctuations. The power index nj^ for the former is given ^] by 

nR = l~6e + 2r] (54) 

where e and rj are defined in Eq. (|40p . We can calculate them based on the approximate solution (jy^^f (t) during the 
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FIG. 6: Time evolution of the scale factor a (t) is shown in logarithmic scale against the cosmic time t. Initial inflation era is 
followed by the almost radiation dominated era. The inset is the later evolution. 



inflation Eq. (P7|) 



3(1 + ^) 



(55) 



which reduces, after use of Eq. ([4T|) . to 



(56) 



where N is the e-folding number defined by Eq. (l4T|) and the next Une. The tensor to scalar ratio r — Pt/Pr — 16e 
becomes 



In our model, they are estimated as 



UR « 0.95 for N = 60, 
r « 4.2 X 10-12 for N = 6G,/ii7inf = 10"'^, k = 0.1. 

These values are well within the observationally allowed range so far [l^ : 



(57) 



(58) 
(59) 



UR = 0.960 ±0.013, 
r < 0.22. 



(60) 
(61) 



We will need to examine other observational tests in the next step in order to confirm the unified BEC cosmology. 
We will save the sequel of this analysis for our future papers. 
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FIG. 7: Time evolution of BEG energy density p^(solid line), boson gas energy density (gray line), and the Hubble parameter 
H (dashed line) in cosmic time t. The original gas density Pg reduces almost completely due to the exponential expansion 
during the inflationary era. In the later stage of the inflationary era, uniform BEG becomes unstable and collapses leaving 
some fraction of p^. At the same time, gradually recovers due to the particle production process caused by the rapid change 
of the background field </> and evetually dominates p^. However pg reduces faster than p^ in the subsequent cosmic evolution. 
Eventually the second instability of BEG takes place. These sequence of BEG and its decay repeats and the effective csmological 
constant reduces toward zero. 

IV. SUMMARY AND DISCUSSIONS 

In the previous papers, we developed the BEC cosmology to describe DE/DM in a unified manner in LA. In this 
paper, we have applied this model to the early universe and unified EA and LA without introducing new ingredients. 

It turns out that the BEC cosmology well fits to describe inflationary dynamics. Being triggered by the condensation 
of bosons in the environment of huge vacuum energy, the inflation naturally initiates. This inflation autonomously 
terminates due to the stagflation stage which inevitably takes place exactly at zero energy density. At the stagflation 
point, particle production and the decay of BEC occur. The former makes the universe connect to the standard hot 
big bang stage and the latter guarantees the vanishingly small cosmological constant after the inflation. Further, we 
have calculated the density fluctuations produced in this model, which turns out to be in the range allowed by the 
present observational data. 

It should be remarked that this stagflation regime is important to solve several problems in the early universe. 
Although this stagflation regime seems to be characteristic to the BEC model, it is not true. The stagflation is 
actually general and robust in various cosmological models provided that the potential of the condensation or the 
inflaton can become negative. It means that there is no necessity to adjust the potential minimum as the zero point of 
energy from the beginning but it is realized autonomously through the stagflation. The stagflation is a normal physical 
mechanism and does not violate any fundamental principles. The total energy density, for example, is guaranteed 
to be non-negative despite the negative potential. Furthermore the potential unbounded from below for a uniform 
mode of the; condensation does not mean the catastrophe of the whole system, but only yields BEC collapse into 
non-uniform modes. 

In order to demonstrate that the stagflation is a general mechanism in the universe, we pick up several typical models 
of inflation and slightly generalize them by adding small negative constant term —V\. The models are the chaotic 
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inflation, the new inflation, and the inflation with exponential potentiaL The potentials of them are generalized to 



Khaotic (0) 
Exponential 



1 



2 J,2 

m (j) 



Vi, 
A 



(62) 

(63) 
(64) 



with TO^ > 0, A > 0, Vi > 0. We numerically calculate the fate of the inflaton fleld (j) {t) in Fig.®. In this figure, the 
columns show, from left to right, inflation models of chaotic, new and exponential potential. The rows show, from 
top to bottom, the potentials V {(j)), time evolutions of BEC order parameter 4>{t), time evolutions of scale factor 
a (t), and time evolutions of energy densities (grey line for Pg (t), solid line for (t)) as well as the Hubble parameter 
H (t) (dashed line). All of the models show stagflation, market by arrows in the figure, without any fine tuning. 
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FIG. 8: Demonstration of the robustness of stagflation in various models of inflation. The columbs show, from left to right, 
chaotic inflation, new inflation and the inflation with exponential potential. The rows show, from top to bottom, the potentials 
V (ip), time evolutions of BEC order parameter <^(t), those of scale factor a{t), and those of energy densities (grey line for 
Pg (t), solid line for p^ (t)) as well as the Hubble parameter H (f)(dashed line). All of them show stagflation, market by arrows, 
without any fine tuning. 



BEC cosmology seems to be able to describe the two accelerating stages of the universe in different ways, one is 
the early inflationary acceleration (EA) and the another is the contemporary mild acceleration (LA). The dynamics 
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of the inflation is triggered by the BEC condensation and the instabihty of BEC around the stagflation point guar- 
antees that the cosmological constant is vanishingly smaU irrespective of the initial vacuum energy. On the other 
hand, the contemporary mild acceleration inevitably takes place from the balance of the condensation speed and the 
potential force. Its attractor nature guarantees the existence of the tiny cosmological constant presently observed. 
EA corresponds to the over- hill regime and LA the mild-inflation regime in the general BEC cosmology [Tj. 

These two acceleration regimes may be connected with each other by multiple collapses of BEC, as suggested in 
Figs.([5][7]). In this paper, we have demonstrated the post-sequence of BEC collapse after inflation and the pre-sequence 
of it before the flnal acceleration. These two sequences may actually be connected to form a single sequence of BEC 
collapse. Unfortunately, the detailed analysis on this possibility and the mechanism of BEC collapse itself go beyond 
the present paper. We will figure out the dynamics of BEC probably by numerical technique and report them in the 
near future. 

Novel properties of BEC cosmology in this paper is realized by the phase transition of BEC, which cannot directly 
described by a unitary Lagrangian dynamics. It will be important to notice that the universe may not be a mechanical 
machine whose time evolution is unitary, but be a whole sequence of phase transitions allowing the emergence of new 
structures in various stages. This paper is also a first proposal for the resolution of the cosmological constant problem 
based on the non-unitary phase transitions including non-uniform modes of condensation. We believe the most 
appropriate method to describe such incoherent evolution is the generalized effective action formalism [ld| on which 
we would like to report comprehensively in our future paper including various applications and verifications of BEC 
cosmology. 
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